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Abstract 



Gravity currents have been studied numerically and experimentally 
both in the laboratory and in the ocean. The question of appropriate 
boundary conditions is still challenging for most complex flows. Grav- 
ity currents make no exception - appropriate, physically and mathe- 
matically sound boundary conditions are yet to be found. This task is 
further complicated by the technical limitations imposed by the current 
oceanographic techniques. 

In this paper, we make a first step toward a better understanding 
of the impact of boundary conditions on gravity currents. Specifically 
we use direct numerical simulations to investigate the effect that the 
popular Neumann, and less popular Dirichlet boundary conditions on 
the bottom continental shelf have on the entrainment and transport of 
gravity currents. 

The finding is that gravity currents under these two different bound- 
ary conditions differ most in the way they transport heat from the top 
towards the bottom. This major difference occurs at medium temper- 
ature ranges. Entrainment and transport at high temperatures also 
show significant differences. 

Key words: Gravity Currents, Entrainment, Transport, Bound- 
ary Conditions, Boussinesq Approximation, Energy Equation 
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1 Introduction 



A gravity or density current is the flow of one fluid within another 
caused by the temperature (density) difference between the fluids P^. 
Gravity currents flowing down a sloping bottom that occur in the 
oceans (e.g., those that flow down the Strait of Gibraltar and the 
Denmark Strait) are an important component of the thermohaline cir- 
culation which in turn is important in climate and weather studies. 

Oceanic gravity currents have been studied in the past, both in 
the laboratory and in the ocean. They have been modelled in various 
ways, starting from the streamtube models, e.g., |3J0] to more recent 
non-hydrostatic Boussinesq equations |SJ [7] • 

The question of appropriate boundary conditions for complex flows, 
such as the gravity currents we consider in this paper, is still challeng- 
ing. This paper is a first step in an effort to bridge the gap between 
observed gravity currents and the assumptions made in modeling them. 
Using ocean observations to develop realistic boundary conditions for 
gravity currents is limited by the available technological means. Thus, 
the current boundary conditions used in the numerical simulation of 
gravity currents are based on physical intuition and mathematical sim- 
plicity and soundness. 

The behavior of gravity currents under two different types of bound- 
ary conditions is studied in this paper. Specifically, the differences be- 
tween gravity currents flowing with Neumann (insulation) and Dirich- 
let (fixed-temperature) bottom boundary conditions are investigated 
via direct numerical simulations. Quite possibly, other boundary con- 
ditions, such as Robin, would be more appropriate. However, given 
the popularity of Neumann [S] IHl El 1X01 CHH an d, to a lesser extent, 
Dirichlet boundary conditions in numerical studies we decided to 
focus first on these two types of boundary conditions. 

Bottom Neumann boundary conditions for temperature may be 
assumed when the material on the continental shelf is a bad thermal- 
conductor (i.e., the current flows over "insulated" rock). Dirichlet 
boundary conditions may be assumed when the material on the conti- 
nental shelf is a good thermal-conductor, making temperature nearly 
constant |12| . 

Dirichlet boundary conditions would be appropriate for the initial 
transient development of gravity currents. For instance, it is known 
that the Red Sea overflow shuts off in the summer ^3] ■ Once this grav- 
ity current starts again, one could expect the temperature difference 
between the bottom layer and the overflow to have a transient impact 
on the mixing near the bottom. Also, such a temperature gradient 
could significantly affect the initial neutral buoyancy level, which is of 
ultimate interest for large-scale ocean and climate studies. Ozgokmcn 
et al. [14] found in numerical simulations that the neutral buoyancy 
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level did not differ significantly from an analytical estimate, which 
does not account for mixing, because the bottom layer properties were 
isolated from vigorous mixing between the gravity current and the am- 
bient fluid, and yet determined the neutral buoyancy level. Thus, any 
mechanism that can affect the properties of near bottom fluid is likely 
to change the neutral buoyancy level, at least during the initial stages 
of the development. This idea will be explored in a future study. 

The rest of the paper is organized as follows: Section|21presents the 
mathematical model used in our numerical study. Section |3| presents 
the numerical model used and the model configuration and parameters. 
Section 01 presents the velocity and temperature boundary and initial 
conditions used in the numerical simulation. Section[S]presents the nu- 
merical investigation of the effect of Neumann and Dirichlet boundary 
conditions on the bottom continental shelf on entrainment and trans- 
port in gravity currents. Five different metrics are used to this end. 
Finally, Sectional presents a summary of our findings and directions of 
future research. 



2 Mathematical Setting 

Consider a two-dimensional gravity current flowing downstream, with 
x as the horizontal (downstream) direction, and z the vertical direction 
(Fig. 1). 

The momentum and continuity equations subject to the Boussinesq 
approximation can be written as: 

a„ A2„. a2„. 

(1) 
(2) 
(3) 

where (x, z) are the two spatial dimensions, t the time, (u, w) the two 
velocity components, p the pressure, g — 9.81 m 2 s~ 1 the gravitational 
acceleration, and i>h and v v the viscosities in the horizontal and vertical 
directions, here assumed to be constants. The way these viscosities 
appear in the simulations is actually via their ratio, r — v v /vh. 

One can argue that, at the horizontal scale being studied (10 km), 
we may assume a constant horizontal viscosity, v^, since its control- 
ling factors, model resolution and the speed of the gravity current (the 
fastest signal in the system) , do not warrant a variable horizontal vis- 
cosity. 

At the scale at which these experiments were conducted, which 
is similar to that in Ozgokmen and Chassignet [H], rotation is not 
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considered important either, i.e., vertical rigidity is not assumed |15| . 
In addition, it is also known that the vertical diffusion in the ocean 
is small |T5]. It is precisely the theory that vertical mixing between 
gravity currents and the ambient fluid is via eddy-induced mixing and 
entrainment, not by diffusion, that led to the assumption that the 
vertical viscosity may likewise be taken to be constant, and assumed 
to be small, in fact. The ratio r — v v /vh is kept at r = 0.5 * 10~ 3 , 
under the range that Ozgokmen and Chassignet estimate the results 
to be insensitive to the values of the vertical viscosity. 
A linear equation of state is used 

p' = -p aT, (4) 

where po is the background water density, a the heat contraction co- 
efficient, and T the temperature deviation from a background value. 
The equation for heat transport is 

dT dT dT d 2 T d 2 T 

— + li— + W— = Kft — + K„ — , (5) 

at ox dz dx z az z 

where k/, and Kjy £1X6 thermal diffusivitics in the horizontal and vertical 
directions, respectively. Nondimcnsionalizing by (x,z) — H(x*,z*), 

(u,w) = (u*,w*), t = f£t*, p = ^p*, T = ATT*, where 

H is the domain depth and AT is the amplitude of the temperature 
range in the system, and dropping (*), the equations become 

du du du dp d 2 u d 2 u 
dt dx dz dx dx 2 dz 2 ' 

dw dw dw dp d 2 w d 2 w 

-w- +u— + w— = -— +RaT+ — T +r—- T , (7) 

at ox oz oz ox 2 - oz 1 

du dw . . 

dx dz ' 

?L u^L dT -p -l ( d2T r d2T \ 

dt dx dz \ dx 2 dz 2 J 

where Ra = (gaATH 3 )/vl is the Rayleigh number (the ratio of the 
strengths of buoyancy and viscous forces), Pr — Vh/i^h the Prandtl 
number (the ratio of viscosity and thermal diffusivity), and r = v v jvy L = 
k v / Kh the ratio of vertical and horizontal diffusivities and viscosities. 

We consider < r < 1. This implies that the viscosity and dif- 
fusivity in cither direction are proportional, and both are small. This 
is in line with the theory that diffusion plays only a small role in the 
vertical mixing between gravity currents and the ambient fluid, and 
the mechanism is via eddy-induced mixing and entrainment. 
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3 Numerical setting 



3.1 Numerical model 

The first consideration taken in passing from the mathematical to the 
numerical model was the small amount of physical dissipation. This 
calls for an accurate representation of the convective operator so that 
the numerical dissipation and dispersion do not overwhelm the assumed 
physical effects. 

Secondly, since small-scale structures are transported with minimal 
physical dissipation, accurate long-time integration is required. Thus, 
high-order methods in space and time are needed. The presence of 
small-scale structures also implies a need for significant spatial resolu- 
tion in supercritical regions, which may be localized in space. 

The numerical method implemented by Nek5000, the software used 
in the simulations, is described below: 

The spatial discretization of ©-© is based on the spectral el- 
ement method 17 , a high-order weighted residual method based on 
compatible velocity and pressure spaces that are free of spurious modes 
Locally, the spectral element mesh is structured, with the solution, 
data and geometry expressed as sums of Nth-order Lagrange polyno- 
mials on tensor-products of Gauss or Gauss-Lobato quadrature points. 
Globally, the mesh is an unstructured array of K deformed hexahcdral 
elements and can include geometrically nonconforming elements. Since 
the solutions being sought are assumed to be smooth, this method is 
expected to give exponential convergence with N, although it has only 
C° continuity. The convection operator exhibits minimal numerical 
dissipation and dispersion, which is required in this study. The spec- 
tral element method has been enhanced through the development of a 
high-order filter that stabilizes the method for convection dominated 
flows but retains spectral accuracy |18|. 

The time advancement is based on second-order semi-implicit op- 
erator splitting methods ^5] GO]- convective term is expressed as 
a material derivative, and the resultant form is discretized via a stable 
backward-difference formula. 

Nek5000 uses an additive overlapping Schwarz method as a pre- 
conditioner, as developed in the spectral element context by Fischer et 
al. |21U22j . This implementation features fast local solvers that exploit 
the tensor-product form, and a parallel coarse-grid solver that scales 
to 1000s of processors |251l2"3] . 

3.2 Model configuration and parameters 

The model domain (see Fig. ^| is configured with a horizontal length 
of L x = 10 km. The depth of the water column ranges from 400 m at 
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x = to H = 1000 to at x = 10 km over a constant slope. Hence the 
slope angle is 9 — 3.5°, which is within the general range of oceanic 
overflows, such as the Red Sea overflow entering the Tadjura Rift [5]. 

Numerical simulations were successfully carried out at a Prandtl 
number Pr = 7 and a Rayleigh number Ra = 320 x 10 6 . 

The numerical experiments were conducted on a Beowulf Linux 
cluster consisting of 9 nodes with 1 Gbps ethernet connectivity. Each 
node has dual Athlon 2 GHz processors with 1024 MB of memory. 
The 2D simulations take approximately 72 hours (simulated to real 
time ratio of ~ 3). 



4 Velocity and temperature boundary and 
initial conditions 



4.1 Temperature boundary conditions 

Inlet (x = 0; 0.6 < z < 1.0): The temperature gradient is one 
of two agents driving the flow, the other being the velocity boundary 
conditions (see Section l4~2"|l . 

The Dirichlet (fixed-temperature) boundary condition is assumed 
at the inlet. This assumption is simple to work with, although its cor- 
rectness ultimately depends on whether or not sea water is an excellent 
thermal conductor |12| . The boundary condition 
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and the initial condition 
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are used (see Fig- EJ| . Figure [3 shows the initial color-coded density 
plot for temperature (red is lowest temperature, blue is highest). 



Bottom (the continental shelf): This is the interface between the 
gravity current and the continental shelf. The actual thermal conduc- 
tivity of the continental shelf will suggest the proper boundary condi- 
tion. Ockendon et al. ^2] prescribe Dirichlet conditions if the boundary 
(the continental shelf) is an excellent thermal conductor; Neumann if 
the region outside of the gravity current (the continental shelf) is a 
very bad thermal conductor; and Robin when the thermal flux at the 
boundary is proportional to the difference between the boundary tem- 
perature and some ambient temperature. This last possibility was not 
studied in this paper, but could be the topic of further study, in the 
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light of the results on incremental transport discussed later in this 
paper. 

Surveys of the literature have not clarified the thermal conductivity 
of the continental shelf. The continental shelf is in the neritic zone, 
so it would contain a mixture of mineral and organic materials. In 
numerical simulations of gravity currents, the most popular choice is 
Neumann boundary conditions QUEUES!- Fewer numerical studies 
employ Dirichlet boundary conditions 1 1 . Other boundary conditions 
(such as Robin) might be more appropriate and should be investigated. 
In this study, we focus on two different sets of boundary conditions: 
Neumann {dT/dn = 0), and Dirichlet (T = 1). 

Dirichlet boundary conditions would be appropriate for the initial 
transient development of gravity currents. For instance, it is known 
that the Red Sea overflow shuts off in the summer JSj- Once this 
gravity current is initiated again, one could expect the temperature 
difference between the bottom layer and the overflow to have a tran- 
sient impact on the mixing near the bottom. Also, such a temperature 
gradient could significantly affect the initial neutral buoyancy level, 
which is of ultimate interest for large-scale ocean and climate studies. 
Ozgokmen et al. |14) found in numerical simulations that the neutral 
buoyancy level did not differ significantly from an analytical estimate, 
which does not account for mixing, because the bottom layer proper- 
ties were isolated from vigorous mixing between the gravity current 
and the ambient fluid, and yet determined the neutral buoyancy level. 
Thus, any mechanism that can affect the properties of near bottom 
fluid is likely to change the neutral buoyancy level, at least during the 
initial stages of the development. This idea will be explored in a future 
study. 

For Dirichlet boundary conditions, a nondimensional temperature 
of one was assumed. 

Numerical experiments were conducted with shelf fixed-temperature 
values < T < 1. T = had the effect of greatly speeding-up the grav- 
ity current, and not allowing the head to form. Thus, the current could 
not entrain ambient water at higher temperature values, contrary to 
experimental and oceanic observations. 

Outlet (x = 10; < z < 1): Note that the outlet is actually a 
sea water to sea water interface, since the region of interest (10 km) is 
only a small portion of the ocean's breadth. The assumed boundary 
condition here is insulation, i.e. dT/dn — 0. (Another possible choice 
is the "do-nothing" boundary conditions |2S|-) Since the outlet is far 
away from the current for most of its journey down the continental 
shelf, the choice of boundary conditions will not influence significantly 
our numerical results. 
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Top (z = 1): This is an air-sea water interface. Presently, insulation 
(dT/dn — 0) is assumed. Since air is a bad thermal conductor, this 
assumption is justified |12| . 

4.2 Velocity boundary conditions 

Inlet [x = 0; 0.6 < z < 1.0): Dirichlct velocity boundary conditions 
are used at the inlet. The velocity profile is given by an "S-shaped" 
polynomial function. The inlet velocity profile is one of the two forces 
that drive the flow, the other being the temperature gradient (see Sec- 
tion 14.11 and Figure OJ) . 

Bottom (the continental shelf): No-slip velocity boundary condi- 
tions (u = 0, w = 0) were used on the bottom. Previous studies 
on gravity currents by Ozgokmen et al. |3] assume no-slip velocity 
boundary conditions on the continental shelf, which is one of the most 
popular assumptions for fluids flowing over a solid boundary. Even in 
their studies on gravity currents over complex topography [TO] , which 
might constitute "roughness" and affect the log layer, Ozgokmen et 
al. still used the no-slip assumption. One should note, however, that 
the question of appropriate boundary conditions flor fluid flows is still 
an active area of research. Other (more appropriate) choices include 
slip-with-friction boundary conditions (23 EH HE! ■ Given the popular- 
ity of Dirichlet boundary conditions in numerical studies of fluid flows, 
however, we decided to focus on them first. 

Interestingly, Haxtel et al. |29| performed experiments on gravity 
currents at high Reynolds numbers for both slip and no-slip, and did 
not observe qualitative changes in the flow structure, although they 
did observe quantitative Reynolds number effects. 

Outlet (x = 10; < z < 1): Free-slip velocity boundary condi- 
tions (dw/dx — 0, dw/dy = 0, du/dx — 0, du/dy — 0) were used 
at the outlet. Another (more appropriate) choice is the "do-nothing" 
boundary conditions |25| . Since the outlet is far away from the current 
for most of its journey down the continental shelf, the choice of veloc- 
ity boundary conditions will not influence significantly our numerical 
results. 

Top (z = 1): Free-slip velocity boundary conditions were used, sim- 
ilar to those used at the outlet. 

Note that there are no boundary conditions imposed on pressure. 
To handle this, Nek5000 uses an operator splitting method, which 
initially uses a "guess for pressure", corrects this guess (but veloc- 
ity is now non-solenoidal) , then corrects the velocity field to impose 
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the divergence-free velocity field. See Rempfer [3U| for a discussion of 
boundary conditions for pressure. 

4.3 Velocity and temperature initial conditions 

The initial velocity field is plotted in Fig. The model is entirely 
driven by the velocity and temperature (Fig.EJ forcing functions at the 
inlet boundary [x — 0), and the temperature gradient. The velocity 
distribution at this boundary matches no-slip at the bottom and free- 
slip at the top using fourth-order polynomials such that the depth 
integrated mass flux across this boundary is zero. 

Other polynomials profiles may be assumed, but care must by taken 
to ensure that 

• the depth-integrated mass flux across the boundary be zero; 

• the assumed temperature profile be consistent with the flow re- 
versal for the velocity boundary condition; and 

• the amplitude of inflow velocity be time-dependent, and scaled 
with the propagation of the gravity current, otherwise there will 
be a recirculation at the inlet (in the case of over-estimation) or 
a thinning-down of the gravity current as it flows down-slope, in 
the case of under-estimation. 

These considerations on initial conditions were taken into account by 
Ozgokmen et al. [HI E3 • 

5 Numerical results 

Gravity currents evolving with two different sets of boundary condi- 
tions (Neumann/insulation and Dirichlet /fixed temperature), are stud- 
ied. Measurements were made of their entrainment, average temper- 
ature, and the heat each transports from the top (warmer water) to- 
wards the bottom (colder water). 

Figure ^ shows the gravity current with Neumann boundary con- 
ditions at 6375 seconds (real time). The head and secondary features 
are very visible at this point. One can see the current entraining some 
of the surrounding fluid via the Kclvin-Hclmholtz instability, which is 
the main mechanism for mixing in a non- rotating fluid [2j • 

FigureElshows the gravity current with Dirichlet conditions, also at 
6375 seconds. One can see a difference in the temperature distributions 
of the two currents, and it is such differences that this study tries to 
quantify. 

The differences between the two currents are most notably in the 
way they transport heat, and in the way they entrain ambient water. 
In order to quantify the differences, five different metrics are used 
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in Sections 15.115.51 These metrics allow a careful assessment of the 
effect of Dirichlet and Neumann boundary conditions on the numerical 
simulation of gravity currents. 



5.1 Entrainment E(t) 

This metric was earlier defined by Morton et al. [SJ. The equivalent 
formulation of Ozgokmen et al. jS] is used here: 

^, * h — h„ , . 

^ = T^£ ■ (10) 

where h is the mean thickness of the current from a fixed reference 
point (which is set at 1.28 km), and h is the mean thickness of the 
tail. The length of the current is taken from the fixed reference point 
to the leading edge of the current. This metric gives us a sense of how 
ambient water is entrained into the gravity current. 

The tail thickness, h Ql is calculated from the flux volumes passing 
through the fixed reference point. The mean thickness, h, is calculated 
from the total volume of the current, from the fixed reference point 
to the nose. The difference in volumes would be accounted for by the 
ambient water entrained. 

The finding here is that the Dirichlet (fixed-temperature) bound- 
ary condition causes greater entrainment than Neumann (insulated) 
boundary condition, by about 14.21%, averaged over the entire travel 
time down the slope. This is actually consistent with the weighted 
temperature metrics (below), which say that the Neumann boundary 
conditions yield a colder current. 

The maximum difference in E(t) between the two currents occurs 
after the head has began to break-down, at around 10,000 seconds 
(actual time). 

5.2 Velocity-weighted temperature 

Similar to the passive scalar calculations of Baringer and Price |32| . 
the velocity-weighted temperature was calculated as 

f (T * u ■ n) dz 
J (u ■ n) dz 

where the integration is from the bottom of the current to its top 
(and the current is understood as having nondimcnsional temperatures 
0.0 < T < 0.8), n is the unit normal (in the x-direction) , and "*" 
denotes multiplication. In this section, and in the remainder of the 
paper, we use nondimensional temperatures. 

The finding is that the Neumann boundary condition yields a cur- 
rent that is 5.61% colder than the Dirichlet case, and the difference 
increases almost linearly with time. 
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5.3 Mean Thickness-weighted Temperature 



Analogous to the velocity- weighted temperature (|llf> . the mean-thickness 
weighted temperature is calculated as 

[ f(T*h)dxdz 

r *<*> - jjhdL ■ ™ 

where the limits of integration are the height (z-direction) and width 
(^-direction) of the current (again, where 0.0 < T < 0.8), and "*" 
denotes multiplication. This is another way of estimating the average 
temperature for the gravity current at any given time. 

Note that metric <|12l) uses the total volume, the same volume calcu- 
lation that is used in calculating h (the mean thickness of the current). 
Hence it is referred to as the mean thickness-weighted temperature. 
On the other hand, the velocity-weighted temperature (fTTf) uses the 
flux volumes instead of the total volume. These flux volumes are what 
is used in calculating h Q (the thickness of the tail). 

By using metric (|12[l . we found that the Neumann boundary con- 
ditions cause the current to be about 9.36% colder than the Dirichlet 
case, slightly more significant than the difference reported by met- 
ric (fTTf) . 

Figure [5] shows the difference in the mean thickness- weighted tem- 
peratures. Just like in Section 15.21 Neumann boundary conditions al- 
ways generate a colder current, and the difference between the average 
temperatures increases linearly with time. 

Both weighted-temperature metrics and (|12|l give results con- 
sistent with the finding for entrainment in Section l5.ll 



5.4 Total heat transport 

Figure [5] shows the initial temperature distribution for both currents. 
The coldest (red) part is at the bottom, and the warmest (blue) part 
at the top. 

The total heat transport for each current was calculated as 

R(t) = J{T*u-n)dz (13) 

at each time-step. R(t) is then normalized with the length of the 
current at that time step. 

The finding is that, on average, the current with Neumann bound- 
ary conditions transports more than Dirichlet case, by 9.34%. 

Figure shows a dramatic rise in the difference in transport be- 
tween the two currents at about 6840 seconds (roughly 1 /3 of the way 
down the domain); the head breaks down at roughly 10,000 seconds 
(real time). 
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5.5 Incremental heat transport 

In the same experiment, repeated calls to the flux-volume calculations 
are done, changing the reference temperature each time. This allows 
the calculation of heat transport for different reference temperatures. 

The data in Section 15 .41 came from the standard reference temper- 
ature of 0.8, which means that only fluid with nondimensional tem- 
peratures less than or equal to 0.8 get into the calculation for R(t) in 
equation (13) above. 

Figure 8 shows the heat transport spectrum plotted for the cur- 
rent with Neumann boundary conditions, at reference temperatures 
0.8,0.7,0.5,0.3 and 0.1. Notice that temperature transport increases 
as the reference temperatures increases. This is because the higher 
the reference temperature, the greater the part of the current that it 
encompasses. Indeed, the standard reference temperature of 0.8 en- 
compasses all temperatures from 0.0 to 0.8. On the other hand, a 
reference temperature of 0.7 encompasses only temperatures from 0.0 
to 0.7. The Dirichlet boundary conditions give a similar spectrum. 

The difference in thermal transport between reference temperatures 
0.5 and 0.7 is calculated, and this gives an estimate of the way temper- 
ature is transported from the top (warmest part) towards the bottom 
(coolest part) of the current. This is equivalent to calculating R(t) 
in equation (13) in the previous section, but only for nondimensional 
temperatures 0.5 < T < 0.7. 

The same calculation is repeated for the middle temperature range 
(from 0.3 to 0.5) and low range (from 0.1 to 0.3). 

The finding here is that, at the middle range, (from 0.3 to 0.5), 
the Dirichlet case transports more by 28.63% (see Fig. |5J). At the 
high range (from 0.5 to 0.7), Dirichlet continues to transport more, by 
11.59%, on average. However, at the low temperature range (from 0.1 
to 0.3), Neumann boundary conditions now give a bigger transport, 
by more than 7.37%. This suggests that there is a major difference in 
the way the two currents transport at medium temperatures, and the 
difference in the incremental transport increases almost linearly with 
time. 

6 Conclusions 

Gravity currents with two types of boundary conditions on the bottom 
continental shelf, Neumann (insulation) and Dirichlet (fixed temper- 
ature), were investigated. The major finding is that the incremental 
temperature transport metric best differentiates the two currents. 

It was noted that the effect is most significant at the medium 
temperature range (0.3 - 0.5), where the difference was 28.63%. In 
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this range, the current with Dirichlet boundary conditions transported 
more heat. 

The other findings (in order of greatest differences) are: The en- 
trainment E{t) metric also differentiates the two currents significantly, 
showing that Dirichlet boundary conditions cause a current to entrain 
by about 14.21% more. The incremental temperature transport at the 
high temperature range showed that the Dirichlet case transports by 
11.59% more. The mean thickness- weighted temperature showed that 
Neumann boundary conditions yield a current that is by about 9.36% 
colder. The transport at standard reference temperature showed that 
the Neumann case transports by 9.34% more. The incremental trans- 
port at the low temperature range showed that the Neumann boundary 
conditions cause about 7.37% greater transport. Finally, the velocity- 
weighted temperature showed that the Neumann case yields a current 
that is by about 5.61% colder. 

Thus, the observed temperature distribution of the gravity cur- 
rent results from a complex interaction of temperature being trans- 
ported from below, and ambient water being entrained from above. 
The present study is just a first step toward a better understanding 
of this complex interaction. In particular, it is found that there can 
be significant differences in the entrainment and transport of gravity 
currents when the most popular Neumann boundary conditions, or 
the less popular Dirichlet boundary conditions are used on the bottom 
continental shelf. 

A better understanding of the impact of both velocity and tem- 
perature boundary conditions on entrainment and transport in gravity 
currents is needed. In particular, slip-with-friction boundary condi- 
tions for the velocity on the bottom continental shelf appear as a more 
appropriate choice. Robin boundary conditions for the temperature 
on the bottom should also be investigated. Different outflow velocity 
and temperature boundary conditions (such as "do-nothing") could 
also yield more physical results, essential in an accurate long time in- 
tegration of the gravity current. All these issues will be investigated 
in futures studies. 
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Figure 1: Set-up 




Figure 2: Color-coded density plot for temperature at time t=0; red is 
coldest, blue is least cold 
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Figure 3: Velocity field at the inlet 
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Figure 5: Gravity current with Dirichlet boundary conditions at time 6375 
seconds 
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Figure 6: Difference in mean thickness-weighted temperature between Neu- 
mann (Neu) and Dirichlet (Dir) cases 



20 



Transport at Ref. Temp = .8 :( Neu - Dir) 




seconds 



Figure 7: Difference in total temperature transport between Neumann and 
Dirichlet cases 
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Figure 8: Transport spectrum for Neumann case for various reference tem- 
peratures 
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Incremental Transport at Medium Temperatures 
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Figure 9: Incremental transport across medium temperature range (0.3 
0.5) 
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